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Surface  Wave  Group  Velocity  Tomography  of  East  Asia 


Introduction 

The  East  Asian  continent  grew  throughout  the  geological  ages  as  a  result  of  successive 
accretions  of  younger  terranes  onto  a  geologically  ancient  (Archean)  core,  the  Siberian  platform 
(e.g.,  Yang  et  al.,  1986).  Figure  1  shows  that  around  the  core  there  are  several  Pre-Cambrian 
geological  terranes,  the  Ordos,  the  Tarim  and  the  Ala  Shan  blocks.  The  Late  Paleozoic  and 
Mid-Mesozoic  South  China  Block  was  most  probably  accreted  to  the  continent  in  Late  Paleozoic 
and  the  Tibetan  Block  was  attached  to  the  continent  in  Mid-  to  Late  Mesozoic.  In  between  these 
older  blocks  were  Paleozoic  collision  zones  (the  Tienshan,  the  Nan-shan  and  the  A-erh-chin 
Shan  [Altyn  Tagh]  Fold  Belts).  The  impingement  of  the  Indian  plate  along  the  Himalayan  front 
that  started  about  fifty  million  years  ago  had  led  to  the  formation  of  Tibet  and  the  ongoing  conti¬ 
nental  tectonics  of  the  whole  East  Asian  area.  While  some  of  the  ancient  blocks,  for  example, 
the  Tarim  and  Ordos,  remain  mostly  passive,  the  fold  belts  were  reactivated  and  became  the 
locus  of  recent  active  tectonics.  The  earthquakes  in  this  area  were  some  of  the  most  energetic 
among  the  continental  events;  the  primary  tectonic  stress  is  thought  to  come  from  the  collision 
south  of  Tibet  (Molnar  and  Tapponier,  1978). 

The  tectonic  units  shown  in  Figure  1  are  drawn  mainly  on  the  basis  of  surface  geology; 
some  of  the  major  units  are  also  clearly  discernible  as  major  topographical  features  (Figure  2). 
But  some  of  the  major  topographical  features  are  not  clearly  expressed  in  the  tectonic  map.  For 
example,  in  Figure  1  Tienshan  is  shown  to  continue  southeastward  of  WMQ  station,  but  Tien¬ 
shan  actually  changes  there  from  high  mountain  ranges  southwest  of  WMQ  to  E-W  oriented 
basin  and  ranges,  with  an  extensive  topographic  low  (Figure  2).  The  lowest  point  is  in  the 
Turfan  basin  (-154  m).  Is  there  a  corresponding  expression  in  the  crustal  and  mantle  structures? 
One  also  wonders  whether  there  is  any  common  deep  structural  expressions  of  the  ancient  blocks 
mentioned  earlier.  Judging  from  the  gradational  changes  of  topography  from  western  China  to 
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eastern  China  (Figure  2),  there  may  also  be  a  corresponding  regional  change  in  crustal  and  upper 
mantle  velocity  structures.  As  far  as  the  Tibetan  plateau  is  concerned,  since  the  manner  in  which 
the  Tibetan  topography  is  supponed  and  the  processes  that  led  to  its  formation  must  be  related  to 
the  deep  stmcture  under  the  plateau,  velocity  structures  under  the  plateau  as  well  as  how  the 
structures  in  the  surrounding  areas  vary  provide  key  infonnation  for  the  contemporaneous  tecton¬ 
ics  of  the  region.  We  conducted  this  surface  wave  tomographic  study  of  Eastern  Asia  to  investi¬ 
gate  the  large  scale  lateral  velocity  variations  in  this  area. 

Reflection  and  refraction  profiles  provide  the  most  detailed  look  of  crustal  structures  of  an 
area.  A  limited  number  of  crustal  reflection  and  refraction  profiles  in  various  parts  of  China  (see 
Wang  and  Mao,  1985,  for  a  summaiy)  have  been  published.  By  far,  the  most  detailed  surve>a 
were  done  in  Nonh  China  Plain,  near  Beijing,  in  conjunction  with  earthquake  studies.  Elsewhere 
there  are  one  or  two  refraction  lines  in  a  large  area,  such  that  the  overall  thickness  of  the  crustal 
thicknesses  are  known.  Several  refraction  and  wide-angle  reflection  studies  were  done  in  south¬ 
ern  Tibet  (Wang  and  Mao,  1985;  Him,  1988).  The  lateral  extent  of  the  plateau,  the  changes  in 
velocity  and  crustal  thickness  within  the  plateau  were  not  mapped  by  such  studies. 

Surface  waves  from  earthquakes,  with  a  wide  period  range,  can  be  used  to  obtain  average 
structures  along  their  travelling  paths  or  tomographic  inversion.  With  the  exception  of  a  two- 
station  dispersion  study  by  Feng  et  al.  (1983),  and  a  tomographic  study  of  Tibet  by  Bouijot  and 
Romanowicz  (1992),  most  of  the  previous  surface  wave  studies  in  East  Asia  were  done  with  data 
external  to  the  region  of  interest.  Feng  et  al.  (1983)  used  data  recorded  on  Kimos  seismographs 
from  stations  within  China  to  derive  Rayleigh  wave  dispersion  in  the  period  range  of  10  to  50 
seconds.  Relatively  few  paths  were  used  in  their  study.  Patton  (1980)  and  Feng  and  Teng 
(1983)  obtained  dispersion  curves  for  various  areas  of  Eurasia  with  Rayleigh  waves  traversing 
through  the  area  by  regionalization;  while  Patton  (1980)  defined  the  sub-regions  based  on  topog¬ 
raphy  and  known  crustal  thickness,  Feng  and  Teng  (1983)  divided  the  region  into  10°  x  10°  grid. 
Velocity  structures  for  the  resolved  regions  were  then  determined  through  inversion.  Chun  and 
Yoshii  (1977)  used  events  on  the  eastern  side  of  the  plateau  and  stations  south  of  the  Himalayas; 
they  aim  their  study  at  Tibet.  Brandon  and  Romanowicz  (1986)  employ  the  "two-event"  tech¬ 
nique  to  determine  dispersion  curves  in  northern  Tibet.  Kozhevnikov  and  Barmin  (1989)  ana- 
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lysed  about  200  records  of  analog  Soviet  stations  and  several  SRO  stations  deployed  in  Asia  to 
obtain  average  Rayleigh  wave  group  velociy  curves  for  several  apriori  determined  tectonic 
regions  of  Eastern  Asia.  These  curves  were  used  by  Kozhevnikov  et  al.  (1992)  to  find  average 
lithosphere  shear  velocity  structure  for  Tibet,  the  mountain  region  of  Southern  Siberia  and  Mon¬ 
golia,  platforms  of  South-eastern  China  and  some  other  regions.  Bouijot  and  Romanowicz 
(1992)  recently  presented  tomographic  images  of  the  Tibetan  Plateau  and  its  vicinity;  they  used 
two  stations  internal  to  China  and  several  SRO  stations  in  Asia  and  also  ANTO  in  Turkey, 

GRFO  in  Germany  and  SSB  in  France. 

For  large  scale  lateral  variations  in  crustal  structures,  the  Bouguer  gravity  map  of  East  Asia 
(USAF,  1971)  provides  a  very  good  view  (Figure  3).  The  most  prominent  anomalies  on  this  map 
are  on  the  Tibetan  plateau,  where  the  anomalies  are  as  low  as  -550  mgal.  The  overall  trend  of 
anomalies  shown  on  the  map  is  the  decrease  in  the  amplitude  of  the  negative  anomalies  toward 
the  east.  More  details  of  this  map  will  be  discussed  when  comparisons  of  the  map  are  made  with 
our  results. 

As  a  result  of  the  establishment  of  high  quality  seismic  stations  in  China  and  its  vicinity, 
the  area  can  now  be  studied  in  more  detail.  The  networks  are  still  too  sparse,  with  station  spac¬ 
ing  on  the  order  of  1000  km.  The  data,  however,  does  provide  an  adequate  basis  for  tomographic 
imaging  of  the  region  as  a  whole.  In  this  paper,  we  present  the  result  of  a  tomographic  study 
using  group  velocities  recorded  at  six  of  the  Chinese  Digital  Seismic  Network  (CDSN)  and  three 
Seismic  Research  Observatory  (SRO)  stations  (see  Figure  2  for  station  distribution).  Our  pri¬ 
mary  purpose  in  this  work  is  to  obtain  an  image,  in  terms  of  group  velocity,  of  the  laterally  heter¬ 
ogeneous  cmstal  and  upper  mantle  structures  of  eastern  Asia.  Group  velocity  is  used  because  we 
want  to  maximize  our  paths  for  better  resolution  and  most  of  the  events  are  too  small  to  have 
dependable  focal  mechanisms. 

We  have  available  also,  through  Prof.  Anatoli  Levshin  a  group  velocity  dataset  determined 
by  Dr.  V,  M.  Kozhevnikov,  using  stations  from  the  former  Soviet  Union  and  China.  By  combin¬ 
ing  this  and  the  CDSN/SRO  dataset  we  can  enlarge  the  area  investigated.  The  preliminary 
results  will  be  shown  in  this  report. 
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The  methodology  used  in  our  tomographic  inversion  was  developed  by  Ditmar  and  Yanov- 
skaya  (1987;  see  also  Levshin  and  Yanovskaya,  1989).  For  each  period,  a  smooth  group  velocity 
distribution  of  the  area  is  covered  by  the  raypaths,  with  its  resolution  (in  km)  depending  on  the 
distribution  of  paths.  The  tomographic  images  of  the  region  as  a  whole  using  Rayleigh  and  Love 
wave  dispersion  curves  show  clearly  the  lateral  variations  in  crustal  and  upper  mantle  structures. 
Although  Tibet  is  by  far  the  most  prominent  feature  in  the  region,  we  are  able  to  resolve  many 
smaUer  features  as  well  in  East  Asia. 

The  tomographic  results  for  the  smaller  dataset  centered  around  China  have  been  sub¬ 
mitted  for  publication  (Wu  and  Levshin,  1993).  Results  for  the  combined  dataset  will  be  sub¬ 
mitted  for  publication,  jointly  with  Dr.  Anatoli  Levshin  after  some  more  tests  and  interpretation 
are  completed.  Initial  results  will  be  reported  at  the  Spring  AGU,  1993  (Levshin,  Wu  and 
Kozhevnikov,  1993). 

Data 

Figure  2  shows  the  locations  of  the  CDSN  and  SRO  stations  and  many  of  the  events  used 
in  this  study.  Because  of  the  wide  dynamic  range  of  the  CDSN  and  the  updated  SRO  seismic 
systems,  although  the  records  stay  on  scale  for  magnitude  7  earthquakes,  surface  waves  from 
=  4.3  can  often  be  used  to  determine  group  velocities  in  the  20-70  seconds  range.  Alto¬ 
gether  100  events,  that  occurred  in  1987,  first  six  months  of  1989, 1990  and  first  half  of  1991, 
were  used;  the  time  span  was  related  to  the  availability  of  data  when  they  were  acquired.  The 
events  chosen  are  located  in  and  around  the  study  area  are  used.  Of  the  500  event-station  paths, 
approximately  360  Love  and  360  Rayleigh  dispersion  curves  were  retained  for  the  final  analyses. 
The  group  velocity  dispersion  curves  were  determined  with  an  interactive  multiple  filter  group 
velocity  program  on  a  workstation,  allowing  rapid  group  velocity  determination  and  visua’  qual¬ 
ity  control.  Dispersion  data  are  discarded  when  the  sonogram  shows  complex  envelop  structures 
along  the  group  arrival.  In  such  cases  we  note  that  the  waveform  is  usually  more  complex  and 
relatively  small;  such  waves  are  probably  radiated  near  the  radiation  pattern  minimum  and  thus 
multipathing  effects  become  pronounced. 
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Because  of  the  relatively  short  path  length  used  in  our  study,  the  effect  if  instrument  group 
delay  is  not  negligible  in  the  determination  of  group  velocities.  McCowan  and  Lacoss  (1978) 
have  shown  that  SRO  instruments  have  delays  ranging  from  about  12  seconds  at  20  seconds  to 
about  24  seconds  at  70  seconds  (without  anti-aliasing  filter).  For  the  CDSN  long  period  channels 
the  group  delays  ranges  from  15  seconds  at  a  period  of  20  seconds  to  35  seconds  at  a  period  of 
30  seconds.  With  the  path  lengths  varying  from  1000  to  4000  km,  neglecting  group  delays  can 
lead  to  an  error  of  several  percent  to  ten  percent.  Figure  4  summarizes  the  group  delays  of  the 
SRO,  CDSN  and  the  modified  SRO  instruments  (MAJ02)  used  in  this  study. 

We  have  recently,  in  cooperation  with  Anatoli  Levshin  of  University  of  Colorado,  joined 
the  dataset  we  have  obtained  using  CDSN  and  SRO  stations  with  analog  stations  from  the  former 
Soviet  Union  and  inverted  for  wider  area  tomography.  The  data  are  from  these  stations  use  a 
mixture  of  instruments,  including  the  intermediate  period  Kimos  seismometers.  The  station 
locations  are  shown  in  figures  displaying  the  results.  The  total  number  of  paths  used  in  this 
combined  study  is  in  excess  of  1200. 

Tomographic  Methodology 

To  inven  surface  wave  group  velocities  we  applied  a  technique  developed  by  Ditmar  and 
Yanovskaya  (1987)  and  Yanovskaya  and  Ditmar  (1990).  This  technique  can  be  considered  a 
generalization  of  the  Backus-Gilbert  inversion  method  (Backus  and  Gilbert,  1968, 1970)  for  2D 
problems  Input  data  for  inversion  are  group  travel  times  r,  for  several  fixed  values 
T„,  m  =  1,  of  period  r  along  given  paths  L,,  y  =  1, and  corresponding  cross¬ 
correlation  matrices  of  travel  time  errors  R,  lr=r„-  Results  of  inversion  are  maps  of  group  veloc¬ 
ity  distribution  C/(0,<}))  I,- ^7-^  and  a  map  of  space  resolution  for  a  given  set  of  paths.  Here  6 

and  ^  are  latitude  and  longitude,  respectively.  The  inversion  procedure  will  be  repeated  for  each 
period  of  interest. 

Let  the  real  distribution  of  group  velocities  be  To  get  a  tomographic  image  of  it 

we  use  a  laterally  homogeneous  initial  model  of  the  area  S  under  study  with  a  constant  group 
velocity  Uq.  The  basic  assumptions  of  the  inversion  technique  are  as  follows: 
i)  Deviations  of  the  real  distribution  of  velocities  from  the  starting  model  are  small,  i.e.. 


(1) 


m,(9,(t))  =  (C/;'(e,<l))-t/-‘)C/o«l 

so  that  we  can  ignore  the  deviations  of  wave  paths  from  great  circles  and,  also,  use  linearized 
inversion  procedures. 

ii)  Taking  into  account  incompleteness  and  inaccuracy  of  the  data  set  we  are  looking  for  a 
smooth  image  w(0,<!))  of  the  real  velocity  perturbations  relative  to  the  starting  model.  To  do  this 
we  introduce  constraints 

J I  Vm(r)  p  dr  -  min  (2) 

s 


and 


^  dm 
y  dn 


)Cs 


=  0. 


(3) 


Here  Q  is  a  contour  surrounding  the  area  S  and  n  is  a  normal  to  Q. 


iii)  The  solution  obeys  constraints  related  to  data 
J  G,{r)mir)dr  =  J  m(s)UQ^ds  =  5r,, 


(4) 

(5) 


Here  C;  is  a  data  kernel  singular  along  the  path  L„  equal  to  zero  outside  the  path  and  obeying  the 


relation 

r 


G.(r)dr  =  to, 

(6) 

\  =  \m{s)Uo'^ds. 

(7) 

Using  the  regularization  technique  (Tikhonov  and  Arsenin,  1976)  in  the  case  of  inaccurate 
data  it  is  possible  to  state  the  problem  of  searching  for  a  solution  as  finding  the  minimum  of  the 
following  functional 
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first  term  in  (8)  is  equal  to  the  total  number  N  of  data.  By  increasing  a  we  impose  stronger 
smoothness  and  decrease  the  resolving  power  of  data. 

The  solution  m(0,(J))  is  found  using  formulas: 


m  =A^bt, 

(9) 

where 

A  =  K\S  +  aJi,r  =  - - to  (5  +  oc;?,)" 

roH5+<xl?,)  Vo 

(10) 

C 

J  Uo 

V 

(11) 

f  C  ds.ds, 

(12) 

Space  resolution  is  determined  using  formula 

R^txpi3J4~A^SA  +2K^A). 

(13) 

Derivation  of  (9-13)  is  given  by  Ditmar  and  Yanovskaya  (1987);  see  also  Levshin  et  al.  (1989). 

The  inversion  proceeds  by  several  steps: 

1.  Transformation  from  spherical  to  Cartesian  coordinates: 

Transformation  from  spherical  coordinates  6,<1)  to  Cartesian 

coordinates  is  done  by  using 

formula  (Yanovskaya,  1982;  Jobert  and  Jobert,  1983): 

X  =/?olntan(9/2). 

(14) 

y 

(15) 

U{x,y)  =  U(B,<^ysme, 


(16) 
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where  /?o  is  the  Earth’s  radius.  This  transformation  does  not  distort  a  velocity  distribution  if  the 

latitude  6  is  not  very  high.  To  make  this  transformation  more  accurate  the  Earth’s  standard  geo¬ 
centric  coordinate  system  is  transformed  first  by  rotation  in  such  a  way  that  the  new  equator 
crosses  the  middle  of  the  territory  under  study  and  the  new  latitude  range  of  wave  paths  is 
narrower  than  the  real  one.  Several  trials  have  demonstrated  that  reasonable  variations  of  the 
new  equator’s  position  do  not  change  results  of  inversion. 

Then  for  a  given  period  T„: 

2.  For  each  T„  the  starting  value  of  a  group  velocity  Uq  is  found  as  an  average  along  all  paths 

S  D, 

J  =  l.J 

Here  D,  is  a  length  of  the  path  L^. 

3.  Functions  f/(0,(J))  and/?(6,(J))  are  found  using  (9-13)  and  transformed  to  initial  coordinate  sys¬ 
tem. 

4.  Steps  2-3  are  repeated  with  different  values  of  the  regularization  parameter.  There  is  an 
option  in  inversion  procedure  taking  into  account  the  presence  of  an  azimuthal  anisotropy  in  the 
Earth’s  model.  The  group  velocity  model  is  constructed  as  a  function  of  coordinates  0,0  and 
azimuth  \j/ 

f^^(e-'}>,V)  =  f/(e.9)(l-fisin2v).  (18) 

The  angle  \j/  and  coefficient  B  are  determined  for  each  point. 

Tomographic  results  for  20°Nto50°N  and  60°Eto  140°E 

The  path  coverage  we  are  able  to  obtain  with  our  present  dataset  is  shown  in  Figure  5.  The 
total  number  of  paths  used  for  each  tomographic  inversion,  the  corresponding  initial  group  velo¬ 
cities  and  the  mean  square  residuals  for  resulting  models  are  presented  in  Table  I.  Figures  6a-e 
and  7a-e  show  the  tomographic  results  for  Rayleigh  and  Love  waves,  respectively,  at  30,  40,  50. 
60  and  70  second  periods.  To  maximize  the  color  contrast  for  each  figure,  we  have  chosen  to  set 
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the  minimum  group  velocity  of  each  figure  to  red  and  the  maximum  to  purple  in  the  rainbow 
color  scale.  The  map  projection  parameters  used  for  these  and  the  topography  map  (Figure  2) 
are  identical,  and  they  can  thus  be  easily  compared.  For  the  tectonic  and  the  Bouguer  anomalies 
the  coordinates  and  the  location  of  the  stations  make  the  comparison  of  these  maps  possible. 
When  viewing  these  maps  it  is  important  to  keep  in  mind  the  following  four  factors.  First,  the 
resolution  maps  shown  in  Figure  6f  and  Figure  6f,  for  Rayleigh  and  Love  waves  respectively;  in 
either  case,  the  resolution  of  our  tomographic  results  is  on  the  order  of  450-700  km  in  the  centrai 
part  of  the  map  and  increases  sharply  toward  the  edge,  where  the  path  coverage  is  poor  (Figure 
5).  Secondly,  since  these  are  maps  of  differences  in  group  velocities,  they  cannot  be  interpreted 
in  terms  of  velocity  structure  differences;  for  example,  low  group  velocities  may  arise  from  a 
combination  of  relatively  low  velocities  in  the  vertical  column  and  a  thick  crust.  Thirdly,  the 
horizontal  wavelength  corresponding  to  30-70  second  waves  are  approximately  95-280  km,  and 
the  image  is  expected  to  be  smoother  for  longer  periods.  Finally,  because  Rayleigh  and  Love 
waves  have  different  eigenfunctions  (Figure  8),  they  sample  the  Earth  differently  even  for  an  iso¬ 
tropic  layered  Eanh;  we  do  not  expect  the  images  for  these  two  w'aves  at  the  same  period  to  be 
identical.  All  results  described  above  were  obtained  assuming  an  isotropic  model  of  the  territoiy' 
under  study.  Inversion  taking  into  '  ccount  the  presence  of  anisotropy  results  in  a  model  w  ith  2^c 
variations  of  velocity  with  azimuth.  The  strike  of  maximum  velocity  direction  strike  is  equal  to 
68'’.  This  model  produces  practically  the  same  group  velocity  maps  and  average  residuals  as  the 
isotropic  one.  We  conclude  that  the  given  set  of  data  does  not  indicate  presence  of  significant 
azimuthal  anisotropy  of  surface  wave  group  velocities  in  the  investigated  regions. 

The  group  velocity  variations  are  remarkably  different  in  different  parts  of  the  study  area. 
To  facilitate  our  description  of  the  maps,  we  shall  divide  this  area  into  three  sub-regions. 

Tibet  and  Southwestern  region 

One  common  feature  clearly  seen  in  all  tomographic  images  (Figures  6  and  7)  are  the 
relatively  low  group  velocities  in  the  western  part  of  the  study  area  when  compared  to  those 
of  the  eastern  part.  The  Tibetan  plateau  becomes  an  outstanding  enclosed  low  group  velocity 
feature  for  Rayleigh  wave  starting  at  40  seconds  (Figure  6b);  at  this  period  it  extends  to  the 
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Pamirs  as  well  as  eastern  Afghanistan  and  nothern  Pakistan.  At  30  seconds  (Figure  6a),  the 
low  group  velocity  area  includes  the  western  Tarim  basin  to  the  northwest  and  to  northeast 
India  and  western  Yunnan.  The  protrusion  toward  station  KMl  is  similar  to  the  shape  of  the 
Yunnan-Tibet  Plateau  in  that  area  (Figure  2).  This  protrusion  is  no  nger  visible  for  periods 
longer  than  40  seconds.  The  area  of  lowest  velocity  is  at  its  maximum  for  40  seconds  and  it 
continues  to  shrink  for  longer  periods,  with  central  Tibet  as  the  core  of  low  group  velocity. 

In  Bourjot  and  Romanowicz’s  (1992)  analysis,  a  similar  velocity  minimum  is  seen  to  persist 
up  to  60  seconds;  in  our  results  we  see  it  clearly  even  at  70  seconds  (Figure  6e),  albeit  the 
area  is  smaller.  Southern  Tibet  and  the  Himalayas  are  areas  of  high  velocity  gradient,  so  are 
the  northern  and  eastern  edge  of  the  Plateau.  The  northwestern  part  of  India,  where  the  reso¬ 
lution  is  good,  appears  as  a  high  group  velocity  region. 

For  Love  wave  tomography  (Figs  7a-e),  Tibet  emerges  as  a  recognizable  low  group 
velocity  feature  for  periods  at  40  seconds.  At  this  period,  the  western  Tarim,  Pamirs  and  the 
Afghanistan-Pakistan  low  that  dominates  the  tomographic  image  at  30  seconds  is  stiU  clear. 
The  Tibet  low,  however,  becomes  the  dominating  feature  at  60  and  70  seconds,  with  the 
overall  shapes  significantly  different  from  that  for  Rayleigh  waves.  But  the  Himalayas  still 
appear  as  a  higher  gradient  zone  and  northwest  India  as  a  relatively  high  velocity  zone. 

Northwestern  region 

Although  the  Tanm  basin  is  a  major  topographic  and  geological  entity,  the  tomo¬ 
graphic  image  does  not  show  it  as  a  distinctive  unit.  Southwest  Tarim  is  a  region  of  low- 
group  velocity  and  forms  a  part  of  the  Tibet  low  group  velocity  anomaly  for  Rayleigh  waves 
at  40,  50  and  60  seconds  (Figures  6b-d).  Otherwise.  Tarim  is  generally  in  the  transition  zone 
from  the  low  group  velocity  region  of  Tibet  to  the  relatively  high  velocity  region  to  the  north. 
Western  Tarim,  however  onsistenUy  shows  up  as  a  low  velocity  unit  in  the  Love  wave 
images  at  all  periods  (Figure  7a-e).  The  Siberian  Shield  in  the  northwestern  comer  of  the 
study  area  is  not  a  well  resolved  region  in  this  study;  Figures  6c-e  do  show  it  as  an  area  of 
relatively  high  velocity. 


As  we  mentioned  earlier,  Tienshan  is  shown  to  continue  from  CIS  Central  Asia  to  the 
northern  edge  of  the  Tarim  Basin  (Figure  1)  on  many  tectonic  maps  (Terman,  1973).  But 
Figures  6a-d  and  7a-e  show  that  the  group  velocities  of  Rayleigh  and  Love  waves  in  the  sec¬ 
tion  of  Tienshan  east  of  LNTS  (see  figures)  is  relatively  high.  Its  possible  significance  will 
be  discussed  later. 

Eastern  region 

Longitude  105°£  is  used  as  the  demarcation  for  the  eastern  and  western  regions  of  this 
study.  A  sharp  group  velocity  transition  exists  near  this  longitude,  especially  in  the  southern 
part  (Figure  6  and  7). 

In  the  eastern  region.  Southeastern  China  is  a  relatively  high  group  velocity  region  at 
all  periods.  The  North  China  Plain  (Figures  1,  6  and  7),  on  the  other  hand,  appears  as  a 
region  with  medium  velocity  at  30  seconds  for  both  Rayleigh  and  Love  waves,  but  becomes  a 
high  velocity  region  at  longer  periods.  Off-shore  the  Ryukyu  backarc  basin  area  has  rela¬ 
tively  low  velocity  at  all  periods  studied,  but  the  Japan  Sea  area  is  shown  as  a  high  velocity 
region  at  30  and  50  seconds  for  Rayleigh  waves  and  a  low  velocity  area  for  other  periods. 

Tomographic  results  for  lO'^NtolO’^N  and  40°Eto  140°E 

In  this  study  we  combine  data  from  the  study  above  and  data  collected  by  Prof.  Anatoli 
Levshin  and  K  from  stations  in  the  former  Soviet  Union.  Although  the  images  in  the  overlap¬ 
ping  area  do  not  differ  too  much  from  those  contained  in  Figures  6  and  7,  the  additional  informa¬ 
tion  in  Figure  8a-e.  We  shall  only  point  out  a  few  key  features. 

Lake  Baikal  and  Mongolia 

Western  and  eastern  Mongolia  is  separated  by  a  ridge  of  high  velocity  that  show  up  in 
all  the  images  shown  in  Figures  8a-e.  It  is  especially  clear  for  T  =  40,  50  and  60.  This  ridge 
coincides  generally  with  the  relative  Bouguer  gravity  high  extending  from  Lake  Baikal  to 
central  Mongolia.  There  is  no  clear  expression  in  the  surface  geology  in  Mongolia,  south  of 
Lake  Baikal  for  this  feature.  As  shown  in  Figures  8b-d  this  high  velocity  ridge  is  surrounded 


on  the  east  and  west  by  two  low  velocity  prongs.  The  nature  of  this  rather  major  feature  is 
most  curious  and  it  could  potentially  have  important  influence  on  reginal  wave  propagation 
through  this  area. 

Siberian  platform 

The  Central  and  the  Western  Siberian  platforms  (Rodriguez,  1969)  are,  in  general,  a 
region  of  high  velocity,  especially  at  longer  periods  (>50  seconds)  as  shown  in  Figures  8c-8e. 
The  relatively  thin  crust  of  the  central  part  of  this  platform  and  the  high  upper  mande  veloc¬ 
ity  under  the  platform  is  probably  responsible  for  the  well  defined  high  velocity  region  at  70 
seconds  (Figure  8  e). 

Discussion 

The  first  order  tomographic  images  presented  in  this  paper  provide  us  synoptic  views  of 
the  lateral  variability  of  the  crust  and  upper  mande  in  East  Asia.  Many  of  the  main  tectonic  units 
(Figure  1)  and  topographic  features  (Figure  2)  can  be  distinguished  quite  well.  Although  we  pre¬ 
fer  not  to  make  further  assumptions,  which  are  necessary  for  the  derivation  of  the  dispersion 
curves  and  then  the  velocity  structures  of  various  regions  from  the  tomographic  results,  our  maps 
do  provide  some  quantitative  measure  of  the  deep  crustal  and  upper  mantle  structures  under  these 
major  features. 

The  Tibet-Pamir  group  velocity  low  dominates  the  tomographic  images  at  periods  greater 
than  40  seconds  (Figures  6b-e).  In  case  of  a  continental  crust,  such  as  that  represented  by  the 
CANSD  model  (Brune  and  Dorman,  1963),  Love  and  Rayleigh  waves  at  40  seconds  or  longer 
have  much  of  their  energy  traveling  in  the  upper  mantle,  and  therefore  sampling  that  part  of  the 
mantle  quite  well.  But  for  a  70  km  Tibetan  crust  (see  Molnar,  1988,  for  wa  summary).  Love  and 
Rayleigh  waves  at  40  seconds  are  trapped  mosdy  in  the  crust  (Figure  9).  These  contrasts  in 
eigenfunctions  are  more  pronounced  at  70  seconds;  at  this  period  the  Tibetan-Pamir  low  is  still 
reflecting  the  thick  low  velocity  crust,  the  relatively  high  velocities  in  eastern  pan  of  China  result 
from  thinner  crust. 
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Judging  from  the  distribution  of  group  velocities  in  the  Tibet-Pamir  area,  the  crustal  thick¬ 
ness  is  most  probably  greatest  in  the  central  pan  of  Tibet.  At  40  seconds,  the  area  of  lowest  Ray¬ 
leigh  group  velocities  includes  the  Kunlun  Terrane,  the  Qiangtang  Terrane  as  well  as  a  part  of 
the  Lhasa  Terrane  nonh  of  the  Himalayas  (Dewey  et  al.,  1988).  As  noticed  by  Bourjot  and 
Romanowicz  (1992),  the  Tibetan  low  velocity  feature  clearly  extends  to  the  southwestern  pan  of 
the  Tarim  Basin  (see  Figures  2  and  6b-c),  where  a  basin  with  thick  (>8  km)  sediments  exists 
(Terman,  1973).  For  Love  waves,  at  40  and  50  seconds  (Figures  7b-c)  the  low  velocity  in  Tibet 
is  not  as  pronounced  as  that  in  Pamir,  but  the  low  velocity  area  is  very  prominent  at  60  and  70 
seconds  (Figures  7d  and  It).  The  low  velocity  area  extends  to  southwest  Tarim  also.  The  low 
velocity  area  in  both  cases  is  surrounded  by  areas  with  high  group  velocity  gradients.  Tarim 
basin  as  a  whole  is  in  the  gradient  zone. 

We  have  noted  earlier  that  the  Tienshan  fold  belt  as  east  of  longitude  87°£  (Figure  1)  is 
noticeably  distinct  from  the  western  part  in  that  whereas  the  western  part  reveals  itself  as  an  area 
with  low  Rayleigh  wave  group  velocity,  the  eastern  pan  is  an  area  of  relatively  high  velocity. 
This  feature  seems  to  be  consistent  with  the  observation  that  the  Bouguer  gravity  low  (Figure  3) 
associated  with  western  Tienshan  (the  -250  mgal  contour)  terminates  there.  Also,  as  noted  ear¬ 
lier,  the  Tienshan  here  is  actually  a  E-W  striking  basin  and  range  province,  with  the  presence  of 
the  sub-sea  level  Turfan  basin  as  the  lowest  point.  Evidently,  this  is  a  deep-seated  feature,  with  a 
thin  crust  underneath,  resulting  perhaps  from  nonh-south  tension.  The  seismicity  of  western  sec¬ 
tion  of  Tienshan  is  rather  high  with  large  thrust  events;  in  contrast  eastern  Tienshan  is  not  very 
seismic. 

The  increase  in  group  velocities  of  both  the  Rayleigh  and  Love  waves  eastward  across 
105°£  is  clear  in  Figures  6b-e  and  7b-e.  The  trend  agrees  generally  with  that  shown  in  the  Bou¬ 
guer  gravity  map  (Figure  3).  In  the  eastern  half  of  the  study  area,  the  relatively  high  velocity 
region  south  and  east  of  Beijing  (the  North  China  Plain)  is  easily  distinguished;  it  is  evidently 
related  to  the  thin  crust  in  that  region,  with  thicknesses  generally  less  than  35  km  (Wang  and 
Mao,  1985);  the  high  group  velocity  here  reflect  the  upper  mantle  structures.  The  North  China 
Plain  is  a  region  of  active  extensional  tectonics  (Nabelek  et.  al.,  1987)  where  many  large  earth- 
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quakes  were  located.  Southeastern  China  is  also  a  region  of  relatively  high  velocities  especially 
at  periods  less  than  70  seconds.  Here  the  thin  crust  (-32  km)  is  probably  the  main  controlling 
factor.  In  contrast  to  North  China  Plain,  this  region  is  not  tectonically  active.  The  Japan  Sea 
area  appears  as  a  high  velocity  region  for  Love  and  Rayleigh  waves  at  40  seconds  (Figures  6b 
and  7b),  but  becomes  an  area  of  relatively  low  velocity  for  longer  period  Rayleigh  waves  (e.g.. 
Figure  6d  and  6e). 

The  interpretation  of  the  new  images  made  with  combined  datasets  is  not  yet  complete. 

One  of  the  obvious  and  very  significant  feature  observed  in  the  images  is  the  high  velocity  ridge 
extending  from  Lake  Baikal  southward  through  Mongolia.  The  nature  of  this  ridge  is  to  be 
investigated  further.  In  view  of  the  fact  that  Lake  Baikal  is  a  young  and  active  tectonic  unit,  this 
southern  extension  may  indicate  a  continuation  of  the  rifting.  Having  started  only  recently,  it  has 
not  yet  had  generated  any  surface  expresion.  Further  studies  are  made  possible  with  the  recent 
establishment  of  the  IRIS/Soviet  network. 

Conclusion 

The  results  of  this  surface  wave  tomographic  inversion  study  provides  clear  images  of  the 
variable  nature  of  the  deep  crustal  and  upper  mantle  structures  under  eastern  Asia  as  a  whole. 

The  surface  wave  tomographic  technique  is  evidently  a  powerful  one  for  this  area  where  a  few 
high  quality  digital  seismic  stations  exist  and  where  there  are  a  large  number  of  moderate  to 
strong  earthquakes  in  and  around  the  study  area.  To  further  refine  the  group  velocity  maps,  it  is 
necessary  to  consider  curved  wave  paths  in  future  work  since  the  velocity  gradients  we  see  in 
these  images  are  quite  high,  with  total  group  velocity  changes  on  the  order  of  ±15%  in  the  study 
area. 
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Table  1. 


Inversion  parameter  and  residuals  for  different  waves  and  periods  for  the  CDSN/SRC  dataset 
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40.1 
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30 

362 
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28.2 

40 
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3.30 

23.7 

50 

362 

3.48 

24.0 

60 

362 

3.61 

21.0 

70 

361 

3.70 

21.9 

Figures 


Figure  1.  Generalized  tectonic  map  of  Eastern  Asia  (after  Terman,  1973). 

Figure  2.  Topography  of  Eastern  Asia  based  on  the  ETOP05  topographic  database.  Locations  of 
stations  and  some  of  the  events  used  in  the  study  are  shown  as  triangles  and  circles  respectively. 
Station  CHTO  is  below  latitude  20  at  longitude  100  E.  Notice  that  many  of  the  tectonic  units 
shown  in  Fig.  1  are  clearly  associated  with  major  topographical  features.  The  Tibetan  Plateau  is 
the  most  prominent  feature  on  this  map.  The  Tarim  Basin  north  of  the  Tibet  where  the  Lop  Nor 
Test  Site  (LNTS)  is  is  bounded  on  the  sorth  by  Tibet  and  the  Tienshan  to  the  north.  Although  the 
eastern  section  of  Tienshan  southeast  of  WMQ  station  is  a  continuation  of  the  western  Tienshan, 
there  is  actually  a  graben  with  its  lowest  point  in  the  Turfan  Basin  (-280m).  The  Szechuan  Basin 
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(red  area  north  of  KMI  station)  is  surrounded  by  1000-2000  mountain  ranges.  The  topography  of 
western  China  is  in  general  much  higher  than  the  eastern  quarter.  The  area  southeast  of  B  J1  sta¬ 
tion  is  the  North  China  Plain. 

Figure  3.  Simplified  Bouguer  gravity  map  of  the  study  area  (USAF, ). 

Figure  4.  Instrument  group  delays  for  World  Wide  Standard  Seismograph  Network  (WWSSN), 
Chinese  Digital  Seismic  Network  (CDSN),  Seismic  Research  Observatory  (SRO)  instruments  at 
TATO  (August  28, 1980  to  mid- 1992)  and  MAJO  (MAJOl  valid  for  the  period  August  23, 1988 
through  August  20,  1990,  and  MAJ02  valid  for  August  20,  1990  through  February  8,  1991). 
Figure  5.  Path  coverage  for  this  study.  Along  most  of  the  paths  both  Love  and  Rayleigh  waves 
are  available.  For  diffemet  periods  the  coverage  varies  slightly. 

Figure  6.  Rayleigh  wave  group  velocity  tomographic  inversion  results  for  (a)  30  seconds,  (b)  40 
seconds,  (c)  50  seconds,  (d)  60  seconds,  and  (e)  70  seconds,  (f)  Resolution  at  50  seconds.  Note 
the  different  scale  for  each  figure. 

Figure  7.  Love  wave  group  velocity  tomographic  inversion  results  for  (a)  30  seconds,  (b)  40  sec¬ 
onds,  (c)  50  seconds,  (d)  60  seconds,  and  (e)  70  seconds,  (f)  Resolution  at  50  seconds.  Note  the 
different  scale  for  each  figure. 

Figure  8.  Rayleigh  wave  group  velocity  tomographic  inversion  results  using  expanded  dataset  for 
(a)  30  seconds,  (b)  40  seconds,  (c)  50  seconds,  (d)  60  seconds,  and  (e)  70  seconds.  Note  the 
different  scale  for  each  figure. 

Figure  9.  Eigenfunctions  for  Rayleigh  (R)  and  Love  (L)  waves  at  40  and  70  seconds  for  two 
extreme  models:  an  approximate  Tibet  model  (T)  with  a  70  km  thick  crust,  and  the  Canadian 
Shield  model  (C)  (Brune  et  al.,  1963).  TL  =  Tibetan  Love  wave,  TR  =  Tibetan  Rayleigh,  CL  = 
Shield  Love  and  CR  =  Canadian  Rayleigh.  For  Rayleigh  waves  only  the  vertical  (radial)  eigen¬ 
functions  are  shown. 
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